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NATIONAL  ADVISORY  COMMITTEE  FOE  AERONAUTICS 
TECHNICAL  NOTE  2693 

A  THEORY  AND  METHOD  FOR  APPLYING  INTERFEROMETRY  TO  THE  MEASUREMENT 
OF  CERTAIN  TWO-DIMENSIONAL  GASEOUS  DENSITY  FIELDS 
By  Walton  L.  Howes  and  Donald  R.  Buchele 


SUMMARY 

A  theory  and  method  are  described  for  the  application  of  interfer¬ 
ometry  to  the  measurement  of  certain  two-dimensional  gaseous  density 
fields.  The  theory  includes  the  effects  of  optical  refraction  upon  the 
observed  interference  pattern.  Exact  equations  denoting  the  relative 
density  difference  corresponding  to  an  observed  interference -fringe 
shift  and  the  optical  distortion  caused  by  refraction  are  derived. 
Corresponding  approximations  in  the  form  of  power  series  expansions  are 
developed  to  permit  practical  application  of  the  theory.  Expressions 
are  derived  which  permit  calculation  of  the  values  of  the  power  series 
coefficients  from  experimental  data. 

The  approximations  were  applied  to  available  interference  data  in 
^  order  to  determine  the  density  distribution  in  a  boundary  layer  formed 

by  supersonic  air  flow  along  a  flat  plate.  Good  agreement  was  obtained 
between  the  density  distribution  calculated  from  the  interference  data 
*  and  that  obtained  from  pressure-probe  measurements. 

Various  methods  for  verifying  the  theory  are  considered  and  an  eval 
uation  process  is  outlined. 


INTRODUCTION 

In  investigations  of  high  subsonic  and  supersonic  gas  flow,  where 
compressibility  of  the  gas  necessitates  the  treatment  of  densi  y  as  a 
variable  various  optical  methods  (reference  l)  have  been  utilized  as  a 

SfS’JST-i®— "one 

in  optical  media  such  as  air  act  in  a  measurable  way  on  light.  More 
over,  optical  methods  of  investigation  possess  the  advantage  of  permitting 
instantaneous  recording  of  the  flow  without  disturbing  the  f low.  Of  the 
optical  methods,  quantitative  investigations  by  the  method  of  interfer 
ometry  have  proved  most  successful.  Interferometric  investigations  of 
flow  fields  have  been  conducted  primarily  with  interferometers  of  the 
Zehnder-Mach  type.  The  principles  of  interferometry  and  the  Zehnder-Mach 
interferometer  have  been  described  elsewhere  (references  2  and  3). 
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The  analysis  by  means  of  interferometry  of  density  distributions 
within  optical  media  is  based  on  the  concept  of  "optical  path  differ¬ 
ences"  of  interfering  light  waves;  the  optical  path  is  defined  as  the 
integrated  product  of  the  physical  path  traversed  by  a  given  wavefront 
and  the  refractive  index  along  that  path.  Thus,  for  convenience,  optical 
path  differences  will  be  expressed  in  terms  of  refractive  index  rather 
than  density.  The  two  quantities  are  related  by  the  Lorentz -Lorenz  law 
(reference  2) .  Also,  for  purposes  of  clarification,  the  mathematical 
development  is  conducted  on  the  basis  of  th'e  ray  theory  rather  than  the 
wave  theory  of  light  propagation.  This  is  possible  because  light  rays 
are  functions  of  the  light  waves  in  that  a  ray  indicates  the  direction 
of  light  propagation,  which,  in  nonpolarizing  media,  is  noimal  to  an 
advancing  wavefront. 


In  cases  such  as  that  considered  herein,  where  the  density  varies  in 
a  direction  inclined  to  the  direction  of  propagation  of  the  light,  the 
mathematical. express ion  of  optical  path  differences  is  complicated  by 
curvature,  or  refraction,  of  the  light  rays, 

Refraction  of  light  in  gas  flow  investigations  by  the  method  of  inter¬ 
ferometry  has  been  considered  by  Weyl  (unpublished),  Bershader  (refer¬ 
ence  4),  Wachtell  (reference  5),  Blue  (reference  6),  and  others.  An 
optimum  wind  tunnel  span  for  minimum  optical  refraction  error  was  deter¬ 
mined  by  F.  J.  Weyl  in  an  unpublished  analysis.  The  first  two  terms 
(constant  gradient  assumption)  of  a  Maclaurin  series  expansion  were  effec- 
tively  utilized  for  refractive  index.  Optical  path  differences  introduced 
outside  the  flow  field  were  neglected.  This  derivation  was  applied  in 
reference  4  to  an  investigation  of  supersonic  air  flow  in  a  channel.  More 
recently,  Weyl's  analysis  was  expanded  by  the  inclusion  of  additional 
optical  path  differences  arising  outside  the  flow  field  and  by  an  expan- 
sion  of  the  Maclaurin  series  representation  of  refractive  index  to  several 
terms  (see  reference  5) .  A  novel  method  of  analysis  was  used  in  refer¬ 
ence  6  to  investigate  boundary  layers  in  supersonic  flow. 


The  present  analysis  attempts  to  include  refraction  effects  in  a 
practical  method  for  evaluating  two-dimensional  gaseous  density  fields 
in  cases  where  the  greatest  space  change  of  density  exists  essentially 
in  one  direction.  The  analysis  was  carried  out  at  the  NACA  Lewis  labor¬ 
atory.  One  coordinate  axis  of  a  Cartesian  coordinate  system  is  oriented 
parallel  to  the  density  gradient.  The  unknown  refractive  index,  or 
density  distribution  is  then  expressed  as  a  Maclaurin  series  with  respect 
STS’  coordinate.  1b.  succeeding  analysis  Is  essentially  an 

and  modification  of  the  analysis  of  reference  5.  On  the  basis  of  the 
succeeding  analysis,  a  procedure  for  evaluating  the  data  obtained  from 

Is  described  and  applied  to  determine  the  density  distri- 
bution  in  a  boundary  layer  formed  by  supersonic  air  flow  along  a 
plate.  Various  methods  are  discussed  for  verifying  the  analysis. 
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PRINCIPLES  OF  METHOD 

A  schematic  diagram  of  the  Zehnder-Mach  interferometer  apparatus  as 
used  in  the  investigation  of  flow-density  fields  is  shown  in  figure  1. 

The  experimental  apparatus  consists  of  a  wind  tunnel  having  a  rectangular 
cross  section  hounded  on  two  sides  by  plane -parallel  glass  windows  and 
a  Zehnder-Mach  interferometer  positioned  with  respect  to  the  wind  tunnel 
so  that  light  (from  a  monochromatic  point  source)  comprising  one  of  the 
divided  collimated  light  beams  is  incident  to  the  wind  tunnel  in  a  direc¬ 
tion  perpendicular  to  the  windows .  The  recombined  light  beams  are  received 
by  a  camera  which  is  focussed  to  image  sharply  some  desired  plane,  for 
example,  plane  O^-Op  :^8ure  within  the  wind  tunnel  and  parallel 

to  the  windows.  When  flow  is  initiated,  the  resulting  nonuniform  density 
field  is  two-dimensional  with  all  density  gradients  essentially  perpen¬ 
dicular  to  the  direction  of  the  incident  light. 

With  the  interferometer  in  correct  adjustment  and  a  condition  of  no 
flow  existing  in  the  wind  tunnel,  an  interference  pattern  consisting  of 
either  a  field  of  uniform  brightness  or  a  series  of  alternate  bright  and 
dark  bands,  or  fringes,  may  be  observed  at  the  image  plane  corresponding 
to  the  plane  I-I  of  figure  1.  The  phase  of  interference  produced  at 
any  point  of  the  image  plane  I-I  by  the  superposition  of  coherent  waves 
is  determined  by  the  optical  path  difference  OPD  of  those  pairs  of  light 
rays  which  originate  at  the  same  point  of  the  light  source  and  intersect 
at  the  same  point  in  the  image  plane. 

When  flow  conditions  are  established  in  the  wind  tunnel,  the  inter¬ 
ference  fringes  observed  at  the  image  plane  may  change  their  position  and 
shape.  The  fringe  shift  at  each  point  of  the  image  is  measured  with 
respect  to  some  reference  pattern,  such  as  the  fringe  pattern  obtained  for 
no-flow  conditions,  or  the  fringe  pattern  representing  some  specific 
region  of  the  flow  field,  for  example,  the  free-stream  region,  for  which 
the  density  can  be  determined.  A  specific  region  of  the  flow  field  is 
adopted  as  the  frame  of  reference  in  the  present  report.  A  detailed 

diagram  of  the  light  paths  is  shown  in  figure  2.  Because  the  observed 

fringe  shifts  are  attributable  to  the  establishment  of  flow  conditions 
in  the  wind  tunnel,  the  light  path  from  the  source  to  the  initial  flow¬ 
boundary  is  the  same  for  both  no-flow  and  flow  conditions .  Therefore, 
the  observed  fringe  shifts  may  be  regarded  as  a  function  of  only  that 
portion  of  the  light  path  which  is  bounded  by  the  initial  flow  boundary, 

B-pB-j-  and  the  image  plane  I-I.  Thus  figure  2  includes  only  the  region 

of  interest  from  the  initial  flow  boundary  Bj-Bj  of  the  wind  tunnel  to 
the  image  plane  I-I.  For  convenience  a  ray  R2,  which  circumvents  the 
wind  tunnel,  is  shown  superimposed  upon  the  wind  tunnel. 


4 


NACA  TN  2693 


Light  Path  in  Plow  Field 

In  order  to  relate  the  unknown  density  field  to  the  observational 
data  secured  by  the  method  of  interferometry,  the  mathematical  relation 
between  the  density  field  and  the  observed  fringe  shifts  must  be  estab¬ 
lished.  Ihis  relation  may  be  obtained  by  treating  the  evaluation  of  the 
two-dimensional  field  as  a  problem  in  one  dimension  whenever  the  density 
gradient  may  be  treated  as  a  function  of  only  one  Cartesian  coordinate. 
Consider  two  Cartesian  coordinate  axes  y  and  r]  both  parallel  to  the 
density  gradient  to  be  investigated  and  lying  in  a  plane  perpendicular  to 
the  direction  of  the  light  path  at  the  initial  flow  boundary  Bj-Bj  . 

(All  symbols  are  defined  in  appendix  A.)  Let  the  y  origin  be  located 
at  a  convenient  point  along  the  density  gradient  and  at  the  initial 
boundary  Bj-Bj  in  figures  1  and  2.  Let  the  q  origin  be  in  the  plane 
Bj-Bj  and  at  that  value  of  y  at  which  any  designated  light  ray 
i  =  0,1,2,.  .  .  enters  the  flow  field.  Let  rj  increase  positively  in 
the  direction  of  increasing  density.  Let  the  units  of  tj  and  y  be 
identical.  Then  tj  and  y  are  related  by  the  transformation 

n  =|y  -  7i|  (!) 

where  y^  refers  to  the  y-ordinate  value  at  which  the  light  ray  i 
enters  the  flow  field,  and 


(2) 


'The  y-coordinate  serves  to  locate  individual  light  rays  with  respect  to 
the  imTmnwn  density  field  whereas  use  of  the  T| -coordinate  permits  the 
mathematical  expression  of  the  ray  paths  in  the  simplest  form. 


In  order  to  interpret  the  unknown  density  distribution  in  terms  of 
the  observed  interference  pattern,  the  paths  of  light  rays  corresponding 
to  interfering  light  waves  must  be  determined .  Hie  light  path  is  given 
by  Fermat* s  principle  as  that  path  for  which  the  time  of  transit  between 
two  points  in  space  attains  a  stationary  value,  usually  a  minimum.  Thus 


/dt  / 


ds 


=  a tat ionary  value 


(3) 


By  definition  the  velocity  v  of  light  propagation  is  related  to  the 
refractive  index  p  of  the  medium  by 

p  =  c/v  (4) 


where  c  is  the  velocity  of  light  propagation  in  vacuum, 
tial  ds  of  physical  path  length  is  given  by 


The  differen- 
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Then  equation  (3)  becomes 


v  =  1/  1  +  <3’X  **  sta"blonary  Talue  (6) 

Let  quantities  possessing  only  the  subscript  i.,  referring  to  a  given 
light  ray,  represent  values  of  the  quantities  at  the  initial  boundary 
Bj-Bj  of  the  unknown  density  field.  Let  quantities  which  are  starred, 
or  subscripted  and  starred,  refer  to  values  of  the  quantities  at  the 
final  boundary  Bji-Bji  of  the  unknown  density  field.  Then  the  solution 

of  the  variation  problem  (equation  (6)),  as  obtained  from  Euler's  equa¬ 
tion  (appendix  B),  is 


which  satisfies  the  boundary  conditions 


0  5  x  1  x* 


0  S  n  *  n* 


and  possesses  a  real  solution  only  when  p(rj)  >  that  is,  when  the 
refractive  index  is  a  monotonic  increasing  function  of  t)  in  which  case 
the  light  ray  curves,  or  refracts,  in  the  xy-plane  toward  the  direction 
of  increased  refractive  index. 

Outside  the  unknown  density  field,  the  refractive  index  is  constant 
within  each  medium;  hence,  it  follows  from  Fermat's  principle  that  the 
light  paths  are  straight  lines.  The  external  path  lengths  can  be  deter¬ 
mined  from  the  geometrical  configuration  of  the  experimental  setup. 

Equation  (7)  may  be  rewritten  as  the  integral  equation 


dx  = 
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the  solution  of  which,  depends  upon  the  function  n .  A  series  solution 
has  been  obtained  from  equation  (7)  by  the  following  method  (reference  5) 
assume  that  the  refractive  index  distribution  throughout  that  region  of 
the  flow  field  traversed  by  any  given  light  ray  may  be  represented  by  the 

Maclaurin  series 


-E 


vhich  converges  to  the  value  a,  "Here  M  ^  refractive  lniei  of  the 

medium  n  la  the  displacement  of  the  light  ray  In  the  wind  tunnel  with 
?  u  _ .  ■ _ _ l _ x. a  +.tae.  t*qtt  ftn+.fiT'pt  the  "hiind  tunne  1 « 


and  by  definition, 


bv 

1 

■“  Vl 

(^») 

\dy  Sy-7l 

Maclaurin  series 


CD 

«  c< 


which  is  assumed  to  converge.  Substitution  of  the  series  and  the 

first  derivative  of  the  series  £  cCTxCT  in  the  differential  equation  (7) 
, .  .  PTnression  in  terms  of  powers  of  x  and  containing  coef- 

ficients  hv  and  c0.  Following  a  eult^de  appr^lmt^n^th^rela- 

tions  of  the  coefficients  bv  and  cCT  of  like  pow 
c0  =  ci  *  c3  =  c5  =  .  .  •  =  c2a-l  -  P 


c2  “  2  bl 


c4al2  blb2 


c6  =  llo  blb22  +  40  bl2^ 


c8  =  5^  blb23  +  slo  *lV5  +  112  \\ 
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where  =  0  because  tj  «  0  at  x  ~  0.  Thus,  the  path  of  a  given 
light  ray  i  in  the  field  \i  is  given  to  a  close  approximation  hy 


X  -  £ 

-  bjX  +  ^2  ^l^*  + 


+  40  ^1 


v6 

X  + 


^5Ck 


^  5  ,  3 

5040  °1°2  +  56o 


*>1  ^2^3  + 


1  v  3, 
112  ^1  ^4 


) 


x8  + 


(13) 


The  location  of  the  light  ray  at  the  final  "boundary  By-B-m  is  found  "by 
replacing  tj  by  «  |y4  -y^|  and  x  by  x  .  Details  of  the  deri¬ 
vation  of  equation  (13)  are  given  in  appendix  B. 

In  addition  to  the  ray  path,  the  angle  a,  through  which  the  light 
ray  is  refracted  in  the  unknown  density  field,  is  also  of  interest.  By 
means  of  the  substitution 

tan  a 


dn 

dx 


(14) 


in  equation  (7)  the  angle  a  may  he  expressed  as 

Hi 

cos  a  -  — 

H 


(15) 


Optical  Path  Relation  of  Interfering  Light  Waves 

The  unknown  density  field  and  the  observed  interference  pattern  are 
related  by  the  optical  path  difference  of  interfering  light  waves.  The 
case  of  a  two-dimensional  density  field  contained  in  a  wind  tunnel  is 
represented  by  figure  2.  In  the  figure,  line  R-^  which  intersects  points 
At,  A2,  A3,  A4,  and  P,  represents  the  path  of  any  arbitrary  light  ray 
which  traverses  the  unknown  density  field.  Superimposed  on  the  diagram 
for  convenience,  the  line  R2  (which  intersects  points  B]_,  B2,  B3 ,  Bg, 

B4,  the  focal  point  F  on  the  lens  L-L,  and  P),  represents  the  path 
of  a  light  ray  which  circumvents  the  unknown  field  and  intersects  the 
light  ray  R^  at  point  P  in  the  selected  image  plane  I-I.  Rays  R]_ 
and  R2  are  constructed  parallel  to  the  optical  axis  at  the  initial 
flow  boundary  Bj-Bj  and  are  assumed  to  originate  at  the  same  point  of 

the  interferometer  light  source.  According  to  the  latter  condition,  the 
light  waves  associated  with  these  rays  are  coherent  and  therefore  may 
interfere . 

The  parallel  construction  of  light  rays  R4  and  R2  implies  that 
the  phase  of  interference  be  constant  over  the  image  plane  I-I,  that  is, 
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that  the  interferometer  by  initially  adjusted  for  obtaining  an  "infinite 
interference  fringe.  However,  it  can  he  easily  shown  that  the  succeeding 
analysis  applies  without  modification  for  any  initial  fringe  adjustment 
and  spacing.  If  an  image  screen  is  located  arbitrarily  at  plane  I_I> 
then  interference  of  the  waves  associated  with  rays  H]_  and  Eg  may  he 
observed  at  point  P.  A  similar  situation  applies  to  all  other  ray  pairs 
emitted  from  the  same  point  in  the  light  source  and  intersecting  at 
plane  I -I. 

According  to  the  geometrical  theory  of  optical  image  formation, 
there  exists  a  conjugate  object  point  associated  with  each  image  point 
such  as  P.  In  the  present  case  it  is  assumed  that  the  optical  layout 
is  such  that  the  object  point  associated  with  P  is  located  in  the 
vicinity  of  the  flow  field.  Kius  ray  B2  intersects,  or  at  least 
appears  to  intersect,  ray  E]_  somewhere  in  the  vicinity  of  the  flow 
field. 


Suppose  that  the  camera  lens  system  L-L  represents  a  perfect 
optical  system  in  the  sense  of  Gaussian  optics.  Then,  in  general,  t  e 
camera  lens  will  produce  a  perfect  point  image  of  a  point  object  when 
the  intervening  object  and  image  spaces  are  represented  by  atmosphere 
under  standard  conditions;  that  is,  all  light  rays  which  intersect  the 
oblect  point  and  traverse  the  lens  will  intersect  at  the  same  point  in 
the  image  space.  If  the  lens  system  is  imperfect  or  if  a  medium,  such 
as  a  glass  wind-tunnel  window,  having  a  refractive  index  which  differs 
from  that  of  standard  atmosphere  is  interposed  between  the  object  and 
image,  or  if  both  conditions  are  present  then  an  imperfect  image  results. 
The  resulting  image  imperfections  are  called  geometrical  aberrations. 

For  the  present,  the  lens  L-L  will  he  regarded  as  perfec  . 

If  the  wind-tunnel  window  is  neglected  and  the  refractive  index  of 
the  object  and  image  spaces  is  considered  to  be  that  of  standard  atmos¬ 
phere,  there  exists  a  conjugate  object  point  which  is  intersects  y 
every* light  ray  which  intersects  P.  This  object  point  is  indicated  by 
the  intersection  of  the  backward  extension  of  ray  Ex  from  point  A3 
with  ray  B2.  The  intersection  is  at  point  PA  lying  in  plane  0A-0A, 

which  represents  the  apparent  Gaussian  object  plane.  Similarly,  all 
other  rays  which  intersect  at  the  image  plane  I-I  also  appear  to  inter¬ 
sect  at  plane  0A-0A . 

How  consider  the  interposition  of  a  wind-tunnel  window  between 
plane  0A-0A  and  the  lens  L-L.  Assuming  that  the  refractive  index  of 
the  surrounding  medium  remains  that  of  standard  atmosphere,  ^ere  ®xis^ 
a  real  Gaussian  object  plane  0R-0R  displaced  from  0A-  A  ^ 
x(t)  independent  of  angle  aa*.  Areal  object  located  in  plane 
(V,-Od  will  in  the  sense  of  Gaussian  optics,  appear  tote  loca^fd  at 
plane  0A-0A  when  the  object  is  observed  through  the  window.  The 

shift  x(t)  is  given  by: 
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x(t)  =  t^l  -  (16) 

where  pa  and  are  the  refractive  indices  of  atmosphere  and  the 

•window,  respectively . 

Next  consider  a  train  of  light  waves  generated  from  point  PA  and 
assume  that  the  refractive  index  throughout  the  object  and  image  spaces 
is  that  of  atmosphere.  According  to  the  theory  of  geometrical  optics, 
the  optical  paths  traversed  by  all  light  waves  from  a  given  object  point 
to  the  conjugate  image  point  of  a  perfect  optical  system  are  equal.  Thus 
the  optical  lengths  defined  by  points  PA,  A3,  A^,  P  and  PA,  B2,  B3, 

B5,  B4,  F,  P  are  equal.  Associated  with  the  light  waves  are  wavefronts; 
a  wavefront  is  defined  as  the  locus  of  points  having  a  constant  optical 
distance  from  PA.  Since  the  refractive  index  is  assumed  to  be  constant 
in  the  object  space,  the  wavefronts  about  PA  are  spheres.  This  con¬ 
dition  follows  from  the  definition  of  optical  path  length.  Consider  the 
wavefront  centered  at  PA  and  intersecting  point  A3.  It  is  apparent 

that  paPAA3  =  paPAB5.  Since  the  optical  path  lengths  from  PA  to  P 
via  points  A4  and  B4  are  also  equal,  it  follows  that  the  optical 
path  length  of  ray  Ex  from  A3  to  P  equals  the  optical  path  length 
of  ray  E2  from  B5  to  P  so  that  these  latter  lengths  can  be 
neglected  because  they  contribute  zero  OPD.  Thus  the  entire  OPD  indi¬ 
cated  by  an  observed  fringe  shift  S  at  point  P  is  represented  by 
the  OPD  of  rays  and  E2  in  the  space  bounded  by  the  initial  flow 

boundary  Bj-Bj  and  the  arc  %%.  The  OPD  of  rays  Ex  and  E2,  as 
indicated  by  a  fringe  shift  Sa  at  point  P,  is  then  given  by 


SaX  " 


b1b2  +  Pw(A2A3  -  t) 


^a  B3B5 


where  Sa,  the  fringe  shift  in  wavelengths  X  of  light,  is  measured 
with  respect  to  the  fringe  located  at  P  when  p  is  replaced  by  pa. 
By  virtue  of  equations  (5)  and  (7),  equation  (17)  may  be  rewritten  as 


SaX  =  W  -  ^a)  x*  +  JI7 


dx  -  +  pw  (A2A3  -  t)  -  pa  B3B5 


where,  in  addition  to  the  substitution  of  equations  (5)  and  (7),  the 
factor  p^x*  has  been  subtracted  from  the  first  term  of  equation  (17), 
added  to  the  second  term,  and  the  terms  have  been  rearranged.  Conse¬ 
quently,  the  first  term  of  equation  (18)  represents  the  OPD  due  purely 
to  the  refractive  index  difference  obtained  for  an  unrefracted  ray 
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■which,  intersects  the  Boundary  B j-Bj  at  A^,  that  is,  at  the  ordinate 
value  Hi  *  0,  or  y  «  yj_.  The  second  and  third  terms  of  equation  (18) 
represent  the  additional  OED  incurred  Between  the  tunnel  Boundaries 

“1  1  _ O _ —  _C» _ _  aJ.  Jt  am  a-C*  v»n  TT  "D 

ll> 


Bt-Bt  and  Bp-B^  as  a  result  of  refraction  of  ray  R-, .  and  the  fourth 


1*  *~jj  -|-  JU  JJ  _ _ 

anil  fifth  terms  indicate  the  OED  introduced  By  the  wind-tunnel  window 
and  the  airspace  Between  the  window  and  camera  lens  B-L,  respectively ,  ^ 
also  as  the  result  of  refraction. 

The  refractive  indices  \x&  of  atmosphere  and  115  of  the  flow  refer¬ 
ence  region,  with  respect  to  which  the  fringe  shifts  are  to  Be  measured, 
are  constant  so  that  a  change  of  refractive  index  within  the  wind  tunnel 
from  that  of  standard  atmosphere  to  that  of  the  flow  reference  region 


yields  a  fringe  shift  Sg  given  By 


S*X  —  ((j.o  -  n„)x 


(19) 


which  corresponds  to  equation  (18)  with  the  exception  that  the  refraction 
terms  are  zero.  Subtraction  of  equation  (19)  from  equation  (18)  gives 

p**  _  _ 

S jX  =  (Hi  -  ^s)x*  +  jir  \  dx  "  x *  +  B3B5 

1  (20) 

where  S-  =  S  -  S§  and  S^  represents  the  observed  fringe  shift  in 

wavelengths  corresponding  to  a  given  light  ray  i  and  measured  with 
respect  to  the  fringe  pattern  representing  the  flow  reference  region. 

Let  the  location  of  the  real  object  plane  Og-Og  Be  determined 
with  respect  to  the  internal  surface  Bji-Bj>  of  the  final  wind  tunnel 
window  By 


X  =  Ex' 


(21) 


Then  it  is  shown  in  appendix  C  By  the  application  of  Snell's  law  of 
optical  refraction  and  consideration  of  the  geometrical  configuration 
in  figure  2  that  the  fringe  shift  S^  in  equation  (20)  is  given  By 
the  exact  expression 


SA  =  (^-(15)2:  + 


\lc  dx  -  M-jX*  +  H-gKx 


^i 


Jo 


4 


a/ 


2  2  *2  , 

^i  ~^i  ~^a 

J  2  ~2  *2 

NV- a  +^i  “^i 


2  2  *2 
"Ri 


r\4^? 


r  1  §  2*2 

A/^a  4Pj  -^j  -(i. 


(22) 
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when  the  camera  lens  system  is  assumed,  to  he  perfect.  Tor  application 
to  practical  problems,  equation  (22)  may  he  simplified  hy  making  suitable 
approximations,  as  shown  in  appendix  C,  with  the  result 


px* 

SiX  =  (Hi  -  H8)x*  +  ~\  H2  cbc  -  Hx*  -  Kx*(p*  -  H)  (23) 


Equation  (23)  is  independent  of  the  wind-tunnel  windows,  their  effect 
upon  the  observed  fringe  shift  S.^  thus  being  negligible.  Assuming 

that  the  series  expansions  -in  equations  (9)  and  (ll)  are  convergent,  it 
is  also  shown  in  appendix  C  that  their  substitution  in  equation  (23) 
yields 


^5 


=  S 


i  x* 


-EE 

v=l  a=  2,4,6 


bv.c, 


OP  00 

vvov  *a  ^  >  > 

x  +  K  Lj  Li 

v=l  cxs:2,4,6 


1 


^VCCTVX 


,  (24) 


where  the  coefficients  cav  may  be  expressed  in  terms  of  the  coeffi- 

cients  ca  by  equating  the  coefficients  of  like  powers  of  x  in  the 
expression 


GO 

£ 

rr=? 


'crv 


xa  = 


(t 


(25) 


Optical  Distortion 


The  refractive  index  difference  -  [15  given  hy  equation  (24) 
is  associated  with  a  given  light  ray  i  which  enters  the  flow  field  at 
the  ordinate  value  y^.  However,  the  fringe  shift  corresponding  to 

the  ray  i  may  appear  to  exist  at  some  ordinate  value  '  y^  +  (A y)^ 
other  than  y.;  that  is,  there  may  exist  a  form  of  optical  distortion 
caused  hy  refraction  of  the  light  rays  in  the  flow  field.  For  example, 
in  figure  2  the  light  rays  E-^  and  Eg  appear  to  originate  at  the 
apparent  object  point  PA,  whereas  the  refractive  index  p^,  indicated 
hy  the  fringe  shift  hy  virtue  of  equation  (24),  exists  at  the 

y-ordinate  value  containing  the  point  A-^  at  which  ray  E-^  enters  the 

flow  field.  Let  the  difference  between  the  y-ordinate  values  contain¬ 
ing  points  PA  and  he  termed  the  "optical  distortion  in  the  object 

space",  denoted  hy  Ay.  Then,  with  respect  to  rays  E^  and  E2,  the 
distortion  may  he  represented  after  consideration  of  figure  2  by 
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(Ar),  =  A3B3  -  (AjBg  +  Bj^) 


(A y)±  -  [kx*  -  x(t)  +  tj  tan  a*  -  £n*  +  y(a{,t)J  (27) 

By  use  of  Snell's  law  and  equations  (15)  and  (16)  it  is  demonstrated  in 
appendix  D  that  equation  (27)  can  be  presented  in  an  exact  form  analogous 
to  that  of  equation  (22)  as  ^ 


(4y)i  -  f01* +  ^  *) 


„*2  -  „,2 


2  2 
^a  +^i 


*  4. 

Tli  +  t 


*2  ,,  2 

H  -  H 


p  9  *; 

V  + 


When  suitable  approximations  are  applied,  equation  (28)  can  be  rewritten 
to  a  close  approximation  as 

(Ay)±  =  Kx*  -  n* 

which  like  equation  (23)  is  independent  of  the  wind-tunnel  windows.  By 
virtue  of  the  differential  equation  (7)  of  the  ray  path,  the  square- 

root  term  of  equation  (29)  may  be  replaced  by  fgj  ,  where  n±  *  1, 

v  /n=ni 


thus  yielding 


(Ay)  =  Kx* 


Then  tj*  is  replaced  by  ^  cCT  xa  and  *>  ^  t]ie  first  derlTa" 

tive  of  Ev'  (assuming  convergence);  equation  (30)  is  now  expressed 


in  its  final  form 


00 

-E 


(Eo  -  1)  c^x 


a=2, 4, 6, . 


Evaluation  Equations 

Equations  (24)  and  (31)  relate  observable  responses  of  individual 
light  rays  to  the  refractive  index  distribution  xn  the  flow  field. 
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Equation  (24)  may  be  readily  transformed  into  an  expression  for  the  rela¬ 
tive  density  change  existing  at  the  value  j±  at  which  a  given  light 
ray  i  enters  the  flow  field  by  applying  the  Glad  stone -Dale  approximation 


(1  a  1  +  kp 


(32) 


of  the  Lorentz -Lorenz  law,  where  p  is  the  gas  density  corresponding  to 
the  refractive  index  p,  and  k  is  a  constant.  Thus,  by  substitution, 
equation  (24)  becomes 


Pj-P5 

P& 


00  00 

EE 


v  — 1  0=2, 4, 6, 


^vcov  *CT 
o+l  X 


+  K 


,4,6, 


bvCov* 


-*o| 


which  may  be  evaluated  if  the  coefficients  bv  are  determined. 


(33) 


Because  of  refraction  of  the  light  rays  traversing  the  flow  field, 
the  observed  interference  pattern  is,  in  general,  a  distorted  represen¬ 
tation  of  the  density  distribution  within  the  field;  the  distortion  of 
the  fringe  pattern  corresponding  to  any  given  light  ray  being  given  by 
equation  (3l) .  In  order  to  determine  the  true  density  field,  the  dis¬ 
tortion  effect  must  be  eliminated  or  at  least  minimized  in  some  manner. 

In  order  to  evaluate  equations  (33)  and  (31),  the  values  of  the 
coefficients  bv,  ca,  and  Cqv  must  be  determined  for  each  ordinate’ 

value  yj_  at  which  any  given  light  ray  i  enters  the  flow  field. 
Actually,  if  the  values  of  the  coefficients  by  are  determined,  the 
corresponding  values  of  the  coefficients  ca  and  Coy  are  given 
automatically  because  the  coefficients  ca  and  Cov  are  functions  of 
bv  according  to  equations  (12)  and  (25),  respectively.  The  coeffi¬ 
cients  bv  may  be  evaluated  if  they  are  expressed  as  functions  of  the 
derivatives  of  the  observed  interference -fringe  shifts  with  respect  to 
y.  For  any  given  light  ray  i,  successive  differentiations  of  the 
assumed  refractive -index-distribution  equation  (9)  yields  equation  (10.), 
where  derivatives  of  the  coefficients  by  are  zero  because  the  coeffi¬ 
cients  by  are  constant  for  any  given  light  ray.  Consider  equation  (24), 
which  characterizes  any  given  light  ray  i,  where  p^  and  bv  possess 

fixed  values.  With  respect  to  the  entire  manifold  of  rays  comprising  the 
light  beam,  p  and  by  (hence  0^)  are  variables.  Thus  equation  (24) 

may  be  rewritten  without  subscripts  i  as 
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—  i-EE  b^~+KEE  ^~( 

X  fc|  4  6,...  Vtsl  a=2,4,6,  . . . 

Differentiation  of  equation  (34)  v  times  with  respect  to  y  yields 


X  +  K 


w  w 

EE 

v  tsl  o  =2 ,4, 6,  •  *  • 


•bvCavx*CT  (34) 


_2l  -iE  f2V  V  ^  x*°  -  kV  V  cavx*CT 

dyWdyV  dyvl4^../+1  44,4,6,..  , 

'  (35) 

Subsequently,  it  follows  from  equation  (10)  that,  at  any  given  ordinate- 
value  yd,  (  .  . 

L  1  X  /  dvs\ 


l3vcav  *a  TJ, 

ZT~  X  ~  h. 

a  +  1 


v=l  a=2,4,6, . . 


*  “v!h**Wv, 


J-7± 


/  oo  ao 

\  v=la=2,4,6, 


hvcCTV 
a  +  1 


vcav  *CT 


oo  oo 

EE byGav  x 

v  —  1  o'— 2 ,4,6,... 


For  simplification,  assume  that  all  terms  on  the  right-hand  side  of 


1 

i — i 

/d”s\ 

bv 

BvTx" 

In  cases  where  the  preceding  assumption  is  invalid,  the  effect  of 
the  neglected  terms  can  he  accounted  for  by  the  application  of  an  iter¬ 
ation  process,  where  the  first  approximation  of  by  is  given  by  equa¬ 
tion  (37) . 

d^S 

It  is  evident  that  the  derivatives  — -  are  determinable  from  the 

dy  . 

observed  interference  pattern.  However,  it  is  important  to  note  that 

aV  O 

the  true  values  of  — rr,  hence  the  desired  values  of  bv,  are  obtained 

dy 

only  when  the  distortion  Ay  is  effectively  zero  throughout  the  region 

of  interest.  However,  it  follows  from  consideration  of  equation  (3l), 

that,  in  general,  no  single  object  plane  of  focus  exists  for  which 
7 

Ay  =  0  for  all  light  rays.  If  the  ray  path  series  13  con~ 

vergent,  then  the  first  nonzero  term  -  for  which  cr  =  2  -  may  generally 
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be  expected  to  be  the  term  of  largest  absolute  magnitude.  In  the  case 
of  equation  (3l)  this  term  may  be  reduced  to  zero  by  selecting  K  =  0.5, 
that  is,  by  selecting  the  object  plane  of  focus  to  be  the  midspan  plane 
of  the  wind  tunnel.  Then  equation  (31)  becomes 


(Ay)1 


£ 

o=4, 6,8, 


o-2  *o 
—  cCTx 


(38) 


Assuming  convergence,  the  magnitude  of  equation  (38)  will  usually  be 
small  enough  so  that  the  shape  of  the  curve  S(y)  will  not  be  changed 

dvS 

appreciably.  Thus  the  values  of  the  derivatives  — y  used  to  compute 

dy 

the  values  of  the  coefficients  bv  will  be  essentially  unchanged  by 
the  residual  distortion.  If  the  curve  S(y)  were  changed  appreciably 
by  consideration  of  equation  (38),  then  an  iteration  process  utilizing 
equations  (37),  (36),  and  (38)  could  be  performed  until  the  desired 
values  of  bv  and.  the  desired  curve  S(y)  were  obtained  such  that  bv 
satisfied  equation  (36)  and  Ay  •*  0  throughout  the  interval  0  =  y  ^  y5 


For  K  =  0.5,  expression  (33)  becomes 


Pj~P5 

P5 


1 

kps 


-£E 

Vssl  <X:=2,4,6,« 


^vcav  *a 

FTT  x 


+ 


00  oo 


"bvcov 


(39) 

Equations  (37),  (12),  (25),  (39),  (38),  and  possibly  (36)  are  suffi 
oient  for  evaluating  an  unknown  gaseous  density  field  satisfying  the 
restrictions  imposed  in  the  derivation. 


Convergence  Considerations 

It  is  impractical  to  investigate  convergence  of  the  series  expan¬ 
sions  and  iteration  processes  analytically  because  the  exact  analytic 

p-pg 

expressions  for  Ay  and  - .  are  unknown.  The  validity  of  the  anal- 

p& 

ysis  can  be  indicated  by  comparing  the  results  obtained  by  the  interfer¬ 
ence  method  with  the  results  obtained  from  other  methods  of  analysis. 

For  example,  the  resultant  distribution  obtained  by  the  interference 
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method  can  he  compared  with  that  given  by  the  theory  for  the  density 
field  or  by  some  other  valid  experimental  method.  Agreement  of  the 
series  solution  with  the  other  methods  of  solution  would  indicate  that 
the  series  are  either  uniformly  convergent  or  asymptotic  representations 
of  the  unknown  quantities.  It  is  possible  that  the  series  might  be 
asymptotic,  although  divergent,  because  the  Maclaurin  series  possesses 
the  form  of  an  asymptotic  series  (reference  7). 

If  the  series  (equations  (38)  and  (39))  are  convergent  or  asymptotic, 
then  in  computations  utilizing  these  series  the  remainder  in  each  case 
is  less  in  absolute  value  than  the  first  series  term  which  is  neglected. 
trhnn  the  number  of  applicable  terms  of  either  series  will  include  all 
successive  terms  for  which  the  coefficients  can  be  calculated  with  a 
reasonable  degree  of  accuracy  up  to,  but  not  including,  the  term  of  mini¬ 
mum  absolute  magnitude. 


EVALUATION  PROCEDURE 


Assume  that  the  interferometer  and  wind  tunnel  are  in  operational 
adjustment  and  that  the  interference  fringes  are  obtainable  for  any 
desired  object  plane  for  which  the  interferometer  camera  may  be  focused 
(before  flow  is  initiated  in  the  wind  tunnel)  by  observing  at  plane  I-I 
(in  fig.  2)  the  image  formed  by  light  reflected  from  an  object  inserted 
into  the  field  of  view  at  the  desired  object  plane  of  focus,  visually  the 
midspan  plane  of  the  wind  tunnel,  for  which  E  =  0.5.  After  the  camera 
is  focussed,  the  object  is  removed,  flow  is  initiated,  and  the  interfero- 
grams  are  obtained  by  means  of  photographic  processes. 

An  interferogram  may  be  evaluated  according  to  any  one  of  several 
schemes  depending  upon  the  initial  choice  of  fringe  orientation  and 
spacing.  The  result  of  the  Interferogram  evaluation  process  is  an 
expression  for  the  fringe  shift  S  as  a  function  of  y  for  the  partic¬ 
ular  experimental  conditions.  These  values  of  .S(y)  can  be  used  In 
the  following  manner  to  determine  the  density  field: 


1.  From  the  experimentally  determined  values  of  S(y),  obtained 
when  K  =s  0.5,  determine  the  values  of  successive  derivatives  gp 


.  d*S  throUghout  the  interval  0  ^  y  ^  y§.  These  values  may  be 
dy2,‘”'dyV 


a-«7  u-j  /  \ 

determined  by  any  one  of  several  methods;  for  example,  S(y)  may  be 
plotted  graphically.  Then  at  any  ordinate  value  y^,  (ay) 


repre- 


-  -  .  ,  y*yt 

sents  the  slope  of  the  curve  S(y).  The  derivatives  can 

.  /  x  ^2o\  y=yi 


plotted.  From  the  curve  values  of 


m 


y=yi 

for  various  values 


y=y. 


* 


* 


* 
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7-J±  can  be  determined  as  before.  The  process  may  be  repeated  to  deter¬ 
mine  successive  derivatives 


dyv 

Alternatively ,  S(y)  may  be  represented  by  a  polynomial  S  =  Pv(y). 
The  values  of  successive  derivatives  f ^ 

Wv 

values  y=yi  may  then  be  calculated  by  taking  successive  derivatives  of 
?v(j)  with  respect  to  y  and  then  ..inserting  the  values  y=yi. 


)y=yi 


at  various  ordinate 


2.  Calculate  values  of  the  coefficients  bv  for  each  light  ray 


from  the  values  of 


w% 


and  equation  (37) . 


3.  Determine  from  the  preceding  results  bv(y)  the  magnitude  of 
the  expression  which  is  contained  in  equation  (36)  but  which  is  neglected 
in  equation  (37).  (The  coefficients  are  given  by  equation  (25)). 

Calculate  the  residual  distortion  given  by  equation  (38),  where  the  coef¬ 
ficients  cCT  are  functions  of  the  coefficients  bv  by  virtue  of  equa¬ 
tions  (12).  If  the  magnitude  of  the  expression  neglected  in  equation  (37) 
iB  of  the  order  of  magnitude  of  the  right-hand  side  of  equation  (37),  or 
if  correction  of  the  residual  distortion  of  the  curve  S(y)  results  in 


an  appreciable  modification  of  the  values  of  the  quantities 
or  if  both  these  conditions  are  true,  then  the  iteration 


fd?8\ 
W  yv/ 


y=yi 

processes  men¬ 


tioned  in  the  section  Evaluation  Equations  may  possibly  be  applied  to 
determine  more  accurate  values  of  (Ay)^  and  bv. 


4.  Calculate  the  relative  density  difference  -  given  by  equa- 

P5 

tion  (39)  for  various  ordinate  values  y=yi. 

P-PS 

5.  Plot  -  against  y,  where  the  observed  values  yn-  are 

PS 

replaced  by  the  distortion-corrected  values  y^  +  (Ay)  ^ .  The  resulting 
curve  represents  the  desired  density  profile. 


NUMERICAL  EXAMPLE 

In  order  to  illustrate  the  application  of  the  equations  and  proce¬ 
dure  and  to  investigate  attendant  conditions,  consider  the  evaluation 
of  the  density  distribution  along  the  density  gradient  within  a  boundary 
layer  due  to  supersonic  air  flow  along  a  flat  plate. 
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The  particular  setup  investigated  (reference  6)  consisted  of  a  flat 
plate  of  3.6-inch  width  contained  in  a  wind  tunnel  of  3. 6 -inch  span 
hounded  by  0.5-inch-thick  glass  windows  having  a  refractive  index  of 
1.52.  The  interferometer  light  source  consisted  of  a  high-pressure 
mercury  arc  in  conjunction  with  an  interference-type  filter  which  trans¬ 
mitted  the  green  spectral  line  of  wavelength  5461  Angstroms  (2.15X10"  in.). 

The  interferometer  camera  was  focussed  to  image  the  midspan  plane  of  the 
wind  tunnel . 

Supersonic  air  flow  having  a  free -stream  Mach  number  of  2.0  at  a 
density  of  6.77X10-4  slugs  per  cubic  foot  was  established.  The  follow¬ 
ing  data  were  obtained  from  an  interf erogram  of  the  boundary  layer 
2.5  inches  downstream  of  the  leading  edge  of  the  flat  plate: 


Fringe  Shift,  S 

Distance  from  plate,  y,  inches 

0 

0.0429 

-.50 

0.0198 

-1.00 

0.0165 

-1.50 

0.0142 

-2.00 

0.0120 

s 

y,  inches 

-2.50 

0.0100 

-3.00 

0.0083 

-3.50 

0.0066 

-4.00 

0.0046 

-4.50 

0.0028 

These  data  are  presented  graphically  in  figure  3.  The  first  two  deri- 
vatives  —  and  SLiL  were  determined  with  a  tangentiometer .  The 

resulting  curves  and-  “ — 3X6  Plotted  in  figures  4  and  5, 

ay  dy 

respectively.  The  quantitative  determination  of  succeeding  derivatives 
was  of  questionable  value.  The  residual  distortion  (indicated  by  equa¬ 
tion  (38)).  the  inherent  error  of  the  original  data  curve ,  and  the  error 
introduced  by  the  measuring  instrument  would  generally  tend  to  introduce 
large  relative  errors  in  the  measured  values  of  succeeding  derivatives. 


— !■  and 

cLy 


The  measured  values  of  ||  and  and  expression  (37)  were 

utilized  to  calculate  the  coefficients1^  and  b2  for  various  values 
of  y  to  obtain  the  curves  b1(y)  and  t2(y),  which  are  also  plotted 
in  figures  4  and  5,  respectively. 

With  respect  to  the  evaluation  equations,  if  refraction  of  the 
light  rays  in  the  flow  field  is  neglected,  then  series  (9)  and  (11) 
refuse  X  p  .  b0  and  r,  -  0,  respectively,  and  the  evaluation  eiua- 

tions  (38)  and  (39)  became,  respectively, 

(Ay),  =  0  (42) 


and 
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pl-p& 

p5 


i  kp6z* 


(43) 


which  represent  first  approximations  of  the  residual  distortion  and 
relative  density  difference.  These  equations  have  been  utilized  by 
several  investigators  with  considerable  success,  especially  when  the 
optical  refraction  was  small.  Equations  (42)  and  (43)  were  used  to 
determine  a  first  approximation  of  the  density  distribution  indicated 
by  the  data  of  the  present  example.  The  resultant  distribution  is 
shown  in  figure  6.  Then,  by  means  of  equations  (37),  (12),  (25),  (38), 

.  dS  d2S 

and  (39)  and  the  derivatives  grr  and  — g,  the  density  distribution 

was  calculated  to  a  third  approximation.  Tbe  appropriate  approximations 
were  determined  as  follows : 

Prom  equation  (37)  it  follows  that 

bn 


X_ /ds\ 

x*  \  dv  / 

\  vy-y i 

(44a) 

1  \  ( d2s\ 

2  X*  l  2  ) 

\dy  /y=yi 

(44b) 

By  virtue  of  equation  (12),  equation  (38)  becomes 

1 


= 


12 


where  only  terms  possessing  coefficients  by,  v  <  2,  hence 

dvS 


(45) 

■a>  o<  4, 


are  included  "because  only  the  values  of  derivatives 


V  9 


2,  were 


determined.  On  this  basis,  the  expansion  of  equation  (39)  is 
pi“p8 


p8 


=  ET  Si  \  ~  |^lc21  +  ^2C22)X*2  +  ft)  ^lc41  +  t2c42)x*4 
K8  L  x 


However,  by  virtue  of  equation  (25) 


°21X*2  +  C41X*^  =  C2X*2  +  C4X*4 


and 
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so  that 


C22X*2  +  °42X*4  •“  + 


°21  =  °2  C41  "  c4  °22  =  0  c42  E  c2 


Therefore,  considering  equation  (12),  the  relative  density  difference 
is  given  hy 

pi"p5  1  1_  ^2  2  .  _L  h?  t*4  ]  / 

57“  =  kpT  (si  X*  “  12  bl  +  30  blb2x  /  ' 


in  the  present  instance. 

The  resultant  density  distribution  (third  approximation)  indicated 
by  equations  (44a),  (44h),  (45),  and  (46)  is  shown  in  figure  6  for  com¬ 
parison  with  the  first  approiimation.  The  second  approximation  is  not 
shown  in  figure  6  because  the  resulting  curve  is  nearly  identical  to 
that  given  by  the  third  approximation..  In  addition,  the  third  approxi¬ 
mation  of  the  density  profile  is  depicted  in  figure  7  for  comparison 
with  the  theoretical  profile  shown  in  reference  6  and  the  experimental 
profile  calculated  in  reference  6  from  total  pressure  measurements. 

The  remainders  of  equations  (45)  and  (46)  were  obtained  hy  expand¬ 
ing  equations  (38)  and  (39),  respectively,  to  include  terms  containing 


with  the  result 


EAy  =  ( blb2  +  20  blb3j  X* 


E  p-p& 


,2.2 

*1*2  + 


Considering  figure  5,  the  maximum  absolute,  value  of  b3  was  found  to  be 
about  5  inches"^.  Thus  the  maximum  absolute  values  of  the  remainders 


were  found  to  be 


<  0.0001  inch  and 


|R  p-ps 


<  0.01. 


In  the  preceding  example,  the  appropriate  evaluation  equations  were 
developed  for  the  case  where  derivatives  of  S  of  order  higher  than  the 
second  are  neglected.  As  noted  previously,  the  calculation  of  values  of 
succeeding  derivatives  would  be  impractical  in  most  cases .  Thus  the 
applicable  formulas  correspond  to  polynomials  possessing  coefficients 
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defined  at  the  initial  "boundary  value  x  *  0,  y  »  y±.  The  evaluation 

*  procedure  corresponds  to  that  which  was  previously  described  with  the 
exception  that  the  former  equations  are  replaced  by  the  expressions 
developed  in  the  preceding  section;  that  is,  the  coefficients  bv  are 
given  by  equations  (44a)  and  (44b),  the  residual  distortion  is  given  by 
equation  (45),  and  the  density  difference  at  each  ordinate  value  yp  is 
calculated  from  equation  (46). 

In  the  preceding  example,  the  principal  basis  for  verification  of 
the  validity  of  the  present  analysis  was  afforded  by  comparison  of  the 
density  field  indicated  by  the  interf erogram  with  that  indicated  by 
theory  and  by  total  pressure  measurements .  The  generally  good  agreement 
of  the  theoretical  and  probe  results  with  the  interferometer  result 
would  tend  to  indicate  that  the  interf erogram  analysis  is  valid  and  prob¬ 
ably  indicates  to  a  good  approximation  the  true  density  distribution  in 
the  present  instance. 

The  maximum  indicated  density  correction  due  to  light  refraction 
amounted  to  a  4-percent  increase  in  the  density  difference  with  respect 
to  free  stream,  or  about  9  percent  of  the  total  density  difference 
across  the  boundary  layer,  as  depicted  in  figure  6.  The  maximum  correc¬ 
tion  occurred  nearest  the  surface  of  the  flat  plate.  This  correction 
resulted  in  considerably  better  over-all  agreement  between  the  density 
field  indicated  by  the  interferogram  and  that  indicated  by  theory  and 
by  pressure  measurements . 

The  maximum  residual  distortion  was  found  from  equation  (45)  to  be 
0.0015  inch  at  a  distance  y  =  0.015  inch  from  the  surface  of  the  flat 
plate.  The  maximum  absolute  remainders  due  to  termination  of  the  series 
were  less  than  0.0001  inch  distortion  and  1-percent  density  change. 

In  the  present  example,  the  previously  mentioned  iteration  processes 
were  neglected  in  the  determination  of  the  coefficients.  bv  and  the  dis¬ 
tortion  Ay.  The  maximum  relative  error  of  the  coefficients  bv 
resulting  from  neglect  of  these  processes  was  estimated  by  means  of  a 
rough  calculation  to  be  about  10  or  20  percent . 

The  present  analysis  assumes  that  the  flow  field  is  perfectly  two 
dimensional  with  the  density  gradient  inclined  essentially  in  one  direc¬ 
tion.  This  condition  is,  however,  only  approximated  in  practice  because, 
for  example,  in  the  illustrative  example  considered  herein,  end  effects 
(reference  4)  serve  to  violate  the  assumed  flow  conditions.  Boundary 
layers  are  formed  along  the  surfaces  of  the  wind-tunnel  windows  and 
adjacent  to  the  windows  and  the  flat  plate.  In  the  example,  the  effec  s 
of  these  deviations  from  the  assumed  conditions  were  assumed  to  be  small. 
However,  preliminary  considerations  indicate  that  "end"  effect  correc- 

*  tions  would  result  in  better  agreement  of  the  profile  indicated  by  inter¬ 
ference  measurements  with  that  indicated  by  theory  and  by  total  pressure 
measurements,  especially  in  the  region  adjacent  to  the  free  stream. 
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Because  no  lens  system  is  absolutely  perfect  it  is  to  be  expected 
that  aberrations  of  the  camera  lens  -will  introduce  additional  optical 
path  differences  and  distortion  (reference  8)  especially  when  the  magni¬ 
tude  of  the  optical  refraction  in  the  flow  field  is  large  and  the  wind 
tunnel -to -camera  lens  distance  is  also  great.  This  effect  is  minimized 
by  the  use  of  a  camera  lens  which  introduces  the  minimum  possible 
geometrical  aberration.  In  the  analysis  and  in  the  example,  the  optical 
path  and  distortion  effects  introduced  by  the  camera  lens  were  regarded 
as  negligible. 


SHIFT  OF  OBJECT  PLANE 

Additional  evidence  of  the  validity  of  the  present  analysis  in  any- 
given  situation  may  be  obtained  by  comparing  the  interference  pattern 
observed  at  an  object  plane  other  than  that  denoted  by  E  =  0.5  with 
that  predicted  by  theory.  The  shift  of  object  plane  introduces  a  double 
effect  upon  the  observed  interference  pattern:  the  phase  of  interference 
anri  the  observed  position  of  the  interference  attributed  to  a  given 
light  ray  are  shifted  because  both  are  functions  of  E  by  virtue  of 
equations  (24)  and  (3l),  respectively.  The  fringe  shift  incurred  in 
translating  the  object  plane  from  the  midspan  plane,  corresponding  to 
E  =  0.5,  to  another  plane  corresponding  to  E  =  E^  may  be  obtained  from 
the  differential  of  equation  (24) ;  that  is, 


dsi  =  “  £  Xy  ]C  ^covX*a+1  a*  (49) 

v*=l  a=2,4,6, . . . 

Similarly ,  the  change  of  distortion  is  obtained  from  the  differential  of 
equation  (31);  that  is, 


* 


dUy^  =  E  a  ca^a  dE  (50) 

0=2,4, 6, ... 

The  shifts  AS1  and  A(Ay)i  due  to  a  shift  AE  of  the  object  plane 

are  determined  by  integrating  equations  (49)  and  (50)  between  the  approp¬ 
riate  limits.  Thus,  for  a  shift  of  object  plane  from  E  =  0.5  to 
E  =  Ex,  AS±  and  A(Ay)i  are  given,  respectively,  by 

oo  oo 

4Si  =  -  X  (%  -  i>  E  E  b,=0V**‘r+1  (51) 

V“1  o^ss2,4,6,  ■  •  • 
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V 


* 


and 

00 

A  (Ay)±  *  (%  -  |)^  <rcax*CT  (52) 

<7=2,4, 6, . . . 

Agreement  of  the  predicted  curve  S(y)  at  E  =  obtained  by  means  of 

equations  (5l)  and  (52)  -with  the  experimental  result  S(y)  would  be 
further  Justification  of  the  analysis  by  means  of  series.  It  is  to  be 
noted  that  in  the  case  where  refraction  of  the  light  rays  in  the  flow 
field  is  neglected,  no  change  in  the  interference  pattern  is  predicted 
when  the  object  plane  is  shifted.  This  condition  is  apparent  from 
equations  (42)  and  (43).  Also,  it  is  to  be  noted  that  the  prediction 
process  is  not  reversible;  that  is,  the  interference  pattern  corres¬ 
ponding  to  the  object  plane  for  which  E  =  0.5  cannot  be  readily  deter¬ 
mined  from  the  interference  pattern  corresponding  to  the  object  plane 
for  which  E  =  E-^.  This  situation  occurs  because  true  values  of  bv 
and  ca  can  be  determined  to  a  close  approximation  only  when  Ay  «  0, 
in  which  case  the  plane  for  which  E  =  0.5  generally  represents  the 
optimum  choice  of  object  plane. 


A  second  set  of  data  S(y)  was  available  for  the  object  plane 
corresponding  to  E  =  1,  that  is,  the  flow -boundary  plane  Bj-Bj  nearest 

the  light  source.  By  means  of  the  values  of  S(y)  for  K=  0.5  and 
equations  (5l)  and  (52),  the  predicted  curve  S(y)  for  E  =  1  was 
calculated,  where,  as  previously,  only  terms  possessing  coefficients 

because  only  the  values  of 

Thus  the  appropriate  expres- 


)v,  v  <  2,  hence  cCT,  a  <  4,  were  included 

V.  W 


dVS 

Lerivatives  — y,  v  <  2,  were  determined, 
3ians  are  ^ 


ASi  =  -  ±  (b2x*3  +  |  bfb2x*5) 


(53) 


and 


A(Ay)i  =  |  bjx*2  +  |  bxb2x*4 


with  the  respective  remainders 


and 


ras  e  -  Tx  ^25bib2  +  To  *1*3) 

(l  *1*2  +  I  *i*s) 


RA(Ay)  20 


(54) 


(55) 


(56) 
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The  maximum  absolute  values  of  the  remainders  were  found  to  he 


<  10 


-1 


and  B 


A  (Ay) 


<10“5  Inch 


The  resulting  curve  S(y)  for  E  s  1  predicted  by  the  present 
theory  is  depicted  in  figure  8  for  comparison  with  the  curve  obtained 
from  experimental  measurements .  Although  the  agreement  of  the  two 
curves  is  not  perfect,  it  is  still  considerably  better  than  that  obtained 
when  refraction  was  neglected,  as  is  evident  from  figure  8.  Apparently 
the  disagreement  between  the  experimental  result  and  that  predicted  by 
the  present  theory  is  caused  principally  by  the  fact  that  the  errors  of 
the  calculated  values  of  the  coefficients  bv  were  relatively  large  when 

the  approximation  (equation  (37))  of  equation  (36)  was  used  alone.  Thus, 
because  the  coefficients  bv,  v  >  0,  comprise  the  major  factors  of  equa¬ 
tions  (51)  and  (52),  the  relative  errors  of  the  quantities  AS^  and 
A(Ay^)  were  also  large.  In  the  case  of  the  evaluation  equation  (39), 
the  effect  of  errors  of  the  calculated  values  of  the  coefficients  bv 
appears  only  in  correction  terms  which  are  small  relative  to  the  primary 
term.  Of  lesser  importance  is  the  fact  that  the  values  of  the  deriva- 


dvS 


were  determined  from  a  curve  S(y)  which  possessed  residual 


tives 

distortion.  Thus,  the  measured  values  of  the  derivatives  would  have  been 
in  error  wherever  the  slope  of  the  experimental  curve  S(y)  deviated  to 
any  great  extent  from  the  slope  of  the  desired  distort ion -free  curve. 


It  is  concluded  that  the  improved  agreement  of  the  predicted  curve 
S(y)  for  K  =  1  (when  refraction  of  the  light  rays  in  the  flow  field 
is  considered)  with  the  experimental  curve  is  further  evidence  of  the 
validity  of  the  series  analysis,  at  least  for  the  present  example. 


GENERAL  DISCUSSION  AND  CONCLUSIONS 

The  present  analysis  represents  an  extension  and  modification  of 
WachteLPs  analysis  with  the  inclusion  of  a  procedure  for  analyzing 
certain  two-dimensional  gaseous  density  fields  by  means  of  optical 
interference . 

The  equation  of  the  light  path  in  the  unknown  density  field  was 
developed  as  a  power  series,  the  coefficients  of  which  are  functions  of 
the  coefficients  of  a  Maclaurin  series  representation  of  refractive 
index.  An  exact  optical  path  difference  equation  was  derived,  which, 
when  re-expressed  in  the  form  of  an  approximation,  permitted  the  deter¬ 
mination  of  the  unknown  density  distribution  from  observed  interference 
fringe  shifts.  A  form  of  optical  distprtion  resulting  from  refraction 
of  the  light  rays  in  the  unknown  density  field  was  considered.  The 


distortion  was  expressed  by  an  exact  equation  and  a  power  series.  The 
power  series  was  utilized  to  determine  an  object  plane  of  focus  for 
which  the  distortion  is  minimized  and  to  calculate  the  residual  distor¬ 
tion.  Differentiation  of  the  optical-path -difference  equation  permitted 
direct  calculation  of  values  of  the  first  few  refractive  index  series 
coefficients . 

In  addition,  a  procedure  for  applying  the  present  theory  was  dis¬ 
cussed  and  applied  to  determine  the  density  distribution  in  a  boundary 
layer  formed  by  supersonic  air  flow  along  a  flat  plate.  The  resultant 
distribution  agreed  reasonably  well  with  that  given  by  theory  and  by 
pressure -probe  measurements.  Further  confirmation  of  the  theory  was 
obtained  by  predicting  relatively  successfully  the  interference  pattern 
produced  for  an  object  plane  other  than  the  optimum  plane. 

It  is  concluded  that  the  present  analysis  will  permit  increased 
accuracy  in  the  determination  of  certain  two-dimensional  gaseous  density 
distributions  by  the  method  of  optical  interference. 


Lewis  Flight  Propulsion  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Cleveland,  Ohio,  January  3,  1952 
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APPENDIX  A 


SYMBOLS 

The  following  symbols  are  used  in  this  report: 


Bp-Bp 


BrBi 


A1,A2,A3,A4  points  intersected  by  light  ray  which  traverses  flow 

field 

Bp-Bp  final  flow  boundary  traversed  by  light  ray 

Bj-Bj  initial  flow  boundary  traversed  by  light  ray 

Bv  binomial  coefficients 

B1,B2,B3,B4,B5  points  intersected  by  light  ray  which  circumvents  flow 

t>v  refractive  index  coefficient 

c  velocity  of  light  in  vacuum 

ca  light  path  coefficient 

F  focal  point  of  camera  lens;  function 

I- 1  image  plane 

K  '  fraction  of  wind  tunnel  span 

k  constant  (0.1166  ft3/slug  for  air  when  X  =  5461A) 

L-L  camera  lens 

0A-0A  apparent  object  plane 

0r-0r  real  object  plane 

OPD  optical  path  difference 

p  image  point 

PA  apparent  object  point 

p(y)  polynomial 

R  remainder 


°A"°A 

0r-0r 
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k1^r2 

S 

s 

t 

V  ■ 

X 

x(t) 


yCoJ^) 

Y±>i\±>  M'i 


„r*  <n*  ..* 

yi^i^i 


a 

A 

T) 

X 

v.or 

P 

Subscripts: 

a 

i 


light  rays 

interference -fringe  shift  (wavelengths  of  light) 
distance  measured  along  path  of  light  ray 
wind-tunnel  window  thickness;  time 
velocity  of  light  in  medium 

distance  from  final  boundary  of  flow  field  to  real 
Gaussian  object  plane 

distance  from  real  Gaussian  object  plane  to  apparent 
Gaussian  object  plane 

coordinates 

! 

displacement  of  light  ray  in  final  wind-tunnel  window 

denotes  value  for  given  light  ray  at  initial  flow¬ 
boundary  traversed  by  ray 

denotes  value  for  given  light  ray  at  final  flow- 
boundary  traversed  by  ray 

angle  measured  with  respect  to  x,  coordinate 

increment  (Ay  »  distortion) 

coordinate  parallel  to  y  coordinate 

wavelength  of  light 

integers 

refractive  index 
density 

atmosphere 

integer,  refers  to  a  given  light  ray 
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w 

8 

v,a 

Superscripts: 


* 


wind-tunnel  window 

flow  reference  region  (free  stream) 

integers 


refers  to  value  at  final  flow-boundary  traversed  by 
light  ray 
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APPENDIX  B 


LIGHT  PATH  IN  PLOW  FIELD 


The  light  path  in  the  flow  field  is  obtained  as  the  solution  of 
the  variation  problem 


Sdt = 


dx  =  stationary  value 


(6) 


where  M-  =  i-i(y).  Equation  (6)  is  of  the  form 

Jat  -Jr(n,sl)  to 

the  solution  of  which  is  given  by  the  Euler  equation 


Sf 

F  -  tj1  -g— p  a  constant 


(Bl) 


(B2) 


where  -n *  =  Thus  the  solution  of  equation  (6)  is 

dx 

1  ■  (tf 

where  \i  a  at  y  =  y^.  Let  the  light  path  be  expressed  by  the 
Maclaurin  series 

00 

CT=0 

where  the  coefficients  ca  are  given  by 


(7) 


(11) 


=  _L  fdan\ 

0  crJ  (  flv-cr  I 

\a*  /  ^=o 


(B3) 


It  follows  that 


oo 

F 


to  =  >  crcaxa-1 
0=1 


(B4) 
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The  assumed  refractive  index  distribution 


-£v 


expression  (B4) ,  and  finally  series  (ll)  are  substituted  in  the  dif¬ 
ferential  equation  (7)  to  obtain  the  expression 


“  (*o) 


bv(cffxff)V  2 


v=l  cr^O 


The  right-hand  side  of  equation  (B5)  may  be  replaced  by  its  binomial 
expansion.  However,  only  the  first  two  terms  of  the  three- term,  expan¬ 
sion  need  be  retained  because  the  coefficient  bg  represents  the 

refractive  index  p-j.  at  tj  =  0.  Since  p  is  very  nearly  unity  over  a 

wide  range  of  temperatures  and  pressures,  the  third  term  of  the  binomial 
expansion  is  negligible.  Therefore,  equation  (B5)  becomes 


cfc„x' 


tv  (C/jX^) 


where  b0  -  1.  .  When  coefficients  of  like  powers  of  x  are  equated, 
the  coefficients  cff  may  be  expressed  as  functions  of  the  refractive 
index  coefficients  bv.  The  first  few  coefficients  cff  are 

c0  “  C1  “  c3  =  c5  = :  ~  c2a_±  =  0 

c2  =  |  bl 


C4  “  12  blb2 

ce  =  m  bi*22  +  jo  bi2 


=8  -  ISO  bi”23  +  m  lv  h;b-5 +  in  bi3  b4 
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Replacing  the  coefficients  ca  in  equation  (ll)  by  the  corresponding 

expressions  in  terms  of  coefficients  bv  results  in  the  following 
expression  for  the  ray  path: 


=  1  blX2 
2  1 


bib2x- 


+  A  b- 


32 


NACA  TN  2693 


APPENDIX  C 

OPTICAL  PATH  RELATION 
Given  the  fringe  shift  expression 

X* 


A 

S±X  =  (pi-p.6)x*  +  -jjy  |i2  dx  -  p±x*  +  (A^-t)  -  M-a  B3B5  (20) 


the  last  two  terms  may  be  expressed  by  consideration  of  figure  2  as 


IV  (A2A3-t)  -  Ayt  ^C0B  ^  -  lj 


(Cl) 


and 


pa  B3B5  =  -  Pa  [kx*  -  X(t)  +  tJ^L  - 
respectively,  where 

X  =  Kx* 


(C2) 


(21) 


The  factors  a*  and  a*  in  equations  (Cl)  and  (C2),  respectively, 
may  be  expressed  in  terms  of  ot*  by  virture  of  Snell’s  law  of  optical 
refraction;  that  is, 

pf  sin  a*  =  p„.  sin  a*  =  pa  sin  a*  (C3) 

and  x(t)  may  be  replaced  by 

-  4  ■  5) 

Then  equations  (Cl)  and  (C2)  become,  respectively. 


r~m 


(16) 


(C4) 


-in2  a* 


Pv  (AgAg-t)  -  p^t 


1 
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* 


Then,  when  (l-cos^  a*)  is  substituted  for  sin^  a*  and  cos  a*  is 
replaced  by  expression  (l5),  that  is. 


v  ^i 

cqs  a*  =  — v 


equations  (C4)  and  (C5)  may  be  rewritten  in  the  form 


respectively.  Substituting  equations  (C6)  and  (C7)  in  equation  (20) 
and  collecting  terms  yield 


The  right-hand  side  of  expression  (22)  consists  of  a  sequence  .of  terms 
comprising  differences  of  quantities  of  nearly  equal  value.  Those  terms 
which  are  included  as  roots  may  be  expanded  by  means  of  the  binomial 
theorem.  Thus 
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where  the  binomial  coefficients  Bv  are  Bq  =  1,  -  t>>  ^2  =  ~  ~q} 

Bj  =  jg,  .  .  .  The  series  (C8)  converge  because  |^2>  Pa2  >>  nf2  “  ^i2* 

Moreover,  because  for  v  >1  the  terms  of  expansions  (C8)  comprise  pro¬ 
ducts  of  differences  of  quantities  of  nearly  equal  value,  such  terms 
may  be  regarded  as  negligible.  Hence  expression  (22)  may  be  reduced  to 

x*  . 

S±X  =  (p-j-p^x*  +  Jf'  li2  dx  -  p±x*  -  |  Kx*(pf2-p12J  (C9) 

where  pa=*  1  except  when  differences  of  quantities  of  nearly  equal  value 
are  concerned.  Furthermore,  because  P-j_” 

p?2  _  ^.2  .  2p±(p£  -H)  (CIO) 

so  that  approximation  (C9)  may  be  rewritten  as 


:  x* 

s±\  =  (p±-p6)x*  +  i  r  p2  dx  -  l^x*  -  Kx*(p*  -M-i)  (23) 

±Jo 

Substituting  expansion  (9)  in  (23)  and  replacing  tj  by  series  (ll)> 
both  expansions  being  assumed  convergent,  yield  the  following  expression 
for  approximation  (23): 
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where  the  third  term  of  the  binomial  expansion  of  p^  is  negligible 
Equation  (Cll)  may  be  rewritten  in  the  form 


where  the  coefficients  cffV  are  expressed  in  terms  of  the  coefficients 
cff  by  equating  the  coefficients  of  like  powers  of  x  in  the 
expression 


Integrating  equation  (C12)  and  solving  for  p^  -  pg  result  in 

OO  OO  CD  OO 

x*a  +  K 


where  the  multiplicative  constant  p-j_  in  the  last  sum  of  -  equation  (C12) 
is  regarded  as  unity. 


36 


NACA  IN  2693 


APPENDIX  D 


OPTICAL  DISTORTION 
From  equation  (27),  that  is, 


(Ay)i  =  [kx*  -  x(t)  +  tj  tan  a*  -  |ji£  +  y(a£,t)J  (27) 

it  follows  by  virtue  of  Snell* s  law  (equation  (C3))  and  equation  (15) 
that 


sin 


<-t.4 


vnf  -  V 


cos  aa  -  J-  V^a2  +  ^i2  "  ^ 


nj2 


tan  a*  =, 

cl 


w.2  2 

Li?  -  U. 


,a2  .♦  ^ 


(Dl) 


Also, 


y(c&t)  -  t  "tan  cty  =  ^  t 


sin  a* 


\  2 

;Y  sin2  aj 


Then,  because  of  equations  (Dl) , 

.  .vT*155 


y« 


I.  2  ,  ..  2  _  .,*2 

+  ^i  ^i 


(D2) 


(D3) 


When  x(t)  is  replaced  "by  expression  (16)  and  equations  (Dl)  and  (D3) 
are  utilized,  equation  (27)  becomes 
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[  M?2-  ^j2 

^a2+^i2"  vf 


~  Pi  +  * 


[  ^2-^j2 

Mv2  +  Hi2  -  ^*2 


Assuming  that  ji±*  p*«  pa=  1,  except  where’  differences  of  quantities 
of  nearly  equal  value  comprise  a  complete  term,  equation  (28)  reduces  to 

(Ay)i  =  Kx*  -  M-±2  -  (29) 

Rather  than  expand  the  square  root  term  of  approximation  (29)  by  the 
binomial  expansion,  which  converges  very  slowly  because  d*  *  |ij_,  the 

term  may  be  replaced  by  (“M  ,  where  1,  which  follows  from 

Vax/Tj=ri* 

equation  (7)  so  that 

Mi  -  **  (§)„.,.  -  li  (30) 

Finally,  because  of  series  (ll)  and  its  derivative,  series  (B4), 


Mi  - 


00 

I 


(K0  -  l)  CaX*a 


CT- —  2  .4.6.. *. 
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Compensator  plates 


BI"|  p  °R 


Image  plane 


Figure  1.  -  Interferometer  apparatus. 
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Figure  5.  -  Variation  of  derivative  of  interf erence-fringe  slope  and  coefficient  b2  with 
distance  from  surface  of  flat  plate.  Fraction  of  wind  tunnel  span,  K,  0.5. 


Density  profiles  indicated  by  interference  measurements 


terference  measurements 
(third  approximation) 
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Figure  7.  -  Comparison  of  density  profiles  indicated  by  various  methods. 
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Figure  8.  -  Predicted  curve  of  fringe  shift  for  fraction  of  wind  tunnel  span  K  =  1  from  data  for  K  =  0.5. 
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